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ABSTRACT 

Several galaxies have exhibited X-ray flares that are consistent with the tidal disrup- 
tion of a star by a central supermassive black hole. In theoretical treatments of this 
process it is usually assumed that the star was initially on a nearly parabolic orbit 
relative to the black hole. Such an assumption leads in the simplest approximation 
to a i -5 / 3 decay of the bolometric luminosity and this is indeed consistent with the 
relatively poorly sampled light curves of such flares. We point out that there is another 
regime in which the decay would be different: if a binary is tidally separated and the 
star that remains close to the hole is eventually tidally disrupted from a moderate 
eccentricity orbit, the decay is slower, typically ~ t~ 12 . As a result, careful sampling 
of the light curves of such flares could distinguish between these processes and yield 
insight into the dynamics of binaries as well as single stars in galactic centres. We 
explore this process using three-body simulations and analytic treatments and dis- 
cuss the consequences for present-day X-ray detections and future gravitational wave 
observations. 

Key words: black hole physics — gravitational waves — hydrodynamics — X-rays: 
general 



1 INTRODUCTION 

In the past few years, several galaxies have exhibited X- 
ray/UV flares consistent with the tidal disruption of a star 
by a supermassive black hole (SMBH; for flare observations 
see |Donley et al.||2002| |Dogiel et al] [20091 |Gezari et aL] 
2009). These candidate disruptions are relevant to the fu- 
eling of some active galactic nuclei (particularly low-mass 
Wang fc Merritt]|2004 1 and contain important in- 



ones; see 



formation about stellar dynamics in the centers of galaxies. 
In addition, they are related to one of the processes be- 
lieved to lead to extreme mass ratio inspirals (EMRIs), in 
which a stellar-mass object spirals into a supermassive black 
hole; EMRIs are thought to be among the most promising 
sources for milliHertz gravitational wave detectors such as 
the Evolved Laser Interferometer Space Antenna (LISA, see 
Amaro-Seoane et al.|2007||2012 l. 



Analyses of stellar tidal disruptions have focused on 
stars whose orbits are nearly parabolic relative to the SMBH 
( [Rees||T988| . In this case, roughly half the stellar material 
becomes unbound and the rest rains down on the SMBH 
with a rate that, for simplified stellar structure, scales with 
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the time t since disruption as M ~ t 5//3 (this is expected 
at late times even for more realistic structure; see |Lodato| 
|et al.|2009| |. 

There is, however, another possible path to disruptions. 
Binaries that get close enough to a SMBH can be tidally sep- 
arated without destroying either star. The result is that one 
star becomes relatively tightly bound to the SMBH whereas 
the other is flung out at high speed. The bound star will un- 
dergo dynamical interactions and its orbit will also shrink 
and circularise due to gravitational radiation. The star may 
eventually be tidally disrupted, but on an orbit that is much 
more bound than in the standard scenario. This will lead to a 



remnant disc of the type analyzed by Cannizzo et al. (1990), 



for which the accretion rate decreases more slowly than in 
the parabolic scenario: M ~ t -1,2 for reasonable opacities. If 
flare light curves are sampled sufficiently these decays could 
in principle be distinguished from each other, which would 
give us new insight into stellar dynamics and the prospects 
for EMRIs. 

Here we present numerical and analytical analyses of 
binary tidal separation and subsequent tidal disruption of 
the remaining star. We note that there exist similar but not 
identic al numerical studies. In particular |Gould fc Quillen] 
|2003| use a mass for the black hole of 3.6 x 10 b M Q but 
show results only for the subset that give captured stars 
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with similar parameters to the observed stars S2-0. Their 
initial binary distributions are similar to ours, although they 
do not examine binaries with initial semi-major axis < 1 



AU and focus on higher masses. Ginsburg & Loeb (2006) 
address a black hole mass of 4 x 10° Mq and their binaries 
are formed of two stars of masses 3M©. They present a 
few sample orbits of captured stars similar to the S-stars, 
but do not give a detailed distribution. |Perets fc Gualandris| 
(2010) also focus on 4 x 10 6 Mq MBHs, and find as expected 
that the captured stars tend to have high eccentricities e > 
0.97, but do not give a periapsis distribution for the stars. 



Madigan et al. ( 2009 \ present in their notable work direct- 



summation A^— body simulations of small discs of stars with 
semi-major axes of 0.026 and 0.26 pc with 4 10 6 Mq MBHs, 
which produced stars with high eccentricities that did not, 
however, enter the region of greatest interest to us. Hence 
we have performed new numerical simulations to explore our 
scenario. 

In § 2 we discuss tidal separations and present our three- 
body simulations of the process. In § 3 we use these results 
as initial conditions and analyze the competition between 
stellar dynamical processes (which can raise or lower the ec- 
centricity) and gravitational radiation (which shrinks and 
circularises the orbit) to determine the mass ranges most 
likely to lead to moderate eccentricities at the point of dis- 
ruption. In § 4 we discuss the tidal process itself, and argue 
that the small but nonzero residual eccentricities mean that 
for sufficiently low-mass SMBHs the star will typically be 
disrupted rather than settling into a phase of steady mass 
accretion onto the SMBH. We present our conclusions in § 5. 



2 BINARY TIDAL SEPARATION 

Tidal separation of binaries by SMBHs was first discussed 
by |Hills| ( |1988[ ) . He suggested that one member of the bi- 
nary would be ejected with a velocity of > 10 3 kms _1 , a 
"hypervelocity star" (HVS); several such objects have now 
been observed (see |Brown et al.|2009| for a discussion of their 
observed properties). The other member would settle into a 
fairly tightly bound orbit around the SMBH; see Mill er et al.| 
( 2005 ) for a discussion in the context of extreme mass ratio 
inspirals into an SMBH. 

To simulate this process we assume a uniform distri- 
bution of pericentre distances between 1 and 700 AU for 
the orbit of the binary around a 10 7 M© MBH. The ini- 
tial orbit is also assumed to be parabolic and to have its 
relative inclination uniformly distributed over a sphere. In 
total 228,000 numerical simulations were conducted using a 



generalized three-body code described by Zare (19741 and 



Aarseth & Zare (19741. This numerical integrator is based 



on Kustaanheimo-Stiefel regularisation of a two-body sys- 



tem, which is described in Kustaanheimo & Stiefel (19651; 
|Aarseth| (p003). The total energy and angular momentum of 
the system are conserved to a high degree of accuracy and 
close encounters between bodies do not induce unphysical 
velocities. 

The resulting distribution for the pericentre distance 
and eccentricity of the captured population, as well as the 
velocity distribution for the star re-ejected into the stellar 
system, are shown in Fig.|l| for an initial internal binary ec- 
centricity of Ei — 0.4 and stellar mass of IMq. To produce 
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Figure 1. Distribution of the pericentre distance and eccentric- 
ity of the captured companion at the tidal separation radius for 
an initial eccentricity of = 0.4 and stellar mass of IMq. Other 
eccentricities do not change significantly the shape of the distribu- 
tion. The red line indicates the maximum pericentre distance for 
which the tidal disruption happens within a Hubble time under 
the influence of gravitational radiation alone. In the limit e — > 1, 
r p approaches the tidal disruption radius, which we display as a 
green line, at at 1.3 AU, although this cannot be seen directly in 
the figure because we are using a resolution of Se = 10 — 4 . 



this figure we chose 10 7 sets of parameters for fixed eccen- 
tricities and drew the semimajor axis of the initial stellar 
binary from a log normal distribution between 0.05 and 10 
AU. This is taken from observations of period distribution 
of binaries in local field stars ( Duquennoy fc Mayor|[T991 1. 



The mean would be about 0.37 AU. 

In the figure we show the resulting probabilities, where 
we plot the probability of finding a captured star with a 
particular pericentre and eccentricity bin given that a binary 
is scattered to within 700 AU of the MBH. The distribution 
of semimajor axes for captured stars is shown in Fig. (|2j) for 
a IMq star that was taken from an initial stellar binary with 
eccentricity ej = 0.0, 0.4, 0.7 or 0.9. 

We now discuss the evolution of the orbits of the stars 
after capture, under the combined influence of two-body re- 
laxation and gravitational radiation. 



3 COMPETITION BETWEEN STELLAR 
DYNAMICS AND GRAVITATIONAL 
RADIATION 

Suppose that a binary has been tidally separated by a close 
passage to a supermassive black hole, but that the remain- 
ing object is outside the tidal radius (i.e., it is not torn apart 
yet). Gravitational radiation will circularise the orbit as it 
shrinks, but dynamical processes can increase the eccentric- 
ity. Eventually, the star will move inside the tidal radius and 
(as we argue in § 4) will probably be tidally disrupted if the 
SMBH is sufficiently low-mass. 

In this section we discuss the dynamics subsequent to a 
tidal separation. We presume that the pericentre of the or- 
bit of the remaining star is outside the tidal radius, so that 
there is no immediate tidal disruption. The star will then 
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Figure 2. Cumulative probability distribution of the semimajor 
axis of the captured star after the tidal separation of the binary 
for a IMq star. The colours denote the initial eccentricity of the 
binary before being disrupted by the MBH, where black (solid 
line) is e; = 0.0, red (dashed line) is = 0.4, green (dot-dashed 
line) is e; = 0.7 and blue (dotted line) is e$ = 0.9. The probabili- 
ties of captures are different in the different eccentricity cases, in 
particular the case e; = 0.9 is easier to capture than the others. 



be subjected to two-body interactions that can change the 
semimajor axis and eccentricity of its orbit. In principle res- 
onant relaxation ( Rauch fc Tremaine|1996 1 could also play a 
role, in particular due to the high eccentricity the orbit has, 
since the component of the torque is linearly proportional 



to eccentricity (Giirkan & Hopman 2007 1, but for the rel- 



evant tight orbits general relativistic pericentre precession 
essentially eliminates this effect (Merritt et al. 20111. We 



will therefore focus exclusively on two-body interactions. 

The two-body energy relaxation time (during which 
the semimajor axis of the orbit will be roughly doubled or 
halved) for a star of mass m moving against a background 
of density p and velocity dispersion a is ( Spitzer|1987 l 



For e ~ 1 this is much shorter than the energy relaxation 
time, hence we will assume that a is fixed throughout. We 
also note that because angular momentum relaxation is a 
random walk process the angular momentum could go up or 
down; if it goes up then nothing interesting happens to the 
star, hence we will consider only the case in which the angu- 
lar momentum and hence the pericenter distance decreases. 

To be more quantitative, let us suppose that we have 
a galactic center with a supermassive black hole of mass 
M with a stellar mass density profile p(r) — po(r/ri n a)~ a 
inside the radius of influence rjnfl = 2GM/ag, where ffo is 
the velocity dispersion in the bulge of the galaxy. The ra- 
dius of influence is by definition the radius inside of which 
the total stellar mass equals the black hole mass, hence 
the normalization is po = . Suppose we make the 

r infl 

simplifying approximation that the velocity dispersion is 
rr(r) = (7o(r/ri n fl) _1 ^ 2 (this scaling is accurate for r <C n n fl 
but not for r ~ ri n fl because of the mass contribution from 
stars). Let us assume in addition an M — op relation of th e 
form M = 10 s M (cr /200 km s" 1 ) 4 ( |Tremaine et al.|2002 |. 



Then r infl w 3M 7 V pc, where M = 10 7 M 7 M Q , and 



t am ~7x 10 11 yr(3 - a^M^mo (r/r infl ) 
(1-e 2 ) 



-3/2 



(3) 



with mo = tti/Mq. f am is the timescale on which two-body 
processes can raise or lower the pericentre distance signifi- 
cantly. Competing against this is the gravitational radiation 
timescale 



£gr ~ 3 x 10 J yr m, 



o 'Mf 2 ( 



1000AU 



(1 



o 2 \ 7 / 2 



(4) 



Over a time t ~ tGR, the orbit shrinks and circularises 
significantly. Setting the two timescales equal to each other 
and noting that the pericentre distance is r p = a(l — e) gives 
a critical pericentre distance of 



0.3 



In A G 2 pm 



(1) 



Here In A ~ 10 is the Coulomb logarithm. For our purposes, 
however, it is not the semimajor axis but the pericentre dis- 
tance that is important, because this is what determines 
whether the star enters the tidal region. It is therefore the 
angular momentum relaxation time that is more relevant. 
For a nearly circular orbit this time is comparable to the 
energy relaxation time, but as we saw in § [2] the initial ec- 
centricity is close to unity in almost all cases. The angular 
momentum of an orbit scales as \J a(l — e 2 ), so an orbit 
with eccentricity e has an angular momentum a factor of 
(1 — e 2 ) 1//2 less than a circular orbit with the same semima- 
jor axis. Two-body relaxation is a diffusive process, hence 
the expected change in energy or angular momentum after 
time t scales as t 1 / 2 . As a result, the angular momentum 
relaxation time is a factor of [(1 — e 2 ) 1//2 ] 2 = 1 — e 2 less than 
the energy relaxation time: 



0.3 



In A G 2 pm 



(1 



(2) 



t « 16 AU (8 x ict 4 )( 2q - 3 )/ 5 



(2a-6)/5 



(5) 



Typical values for r P: CI u can be read directly off of the sim- 
ulations. For one that can decay faster than a Hubble time 
it is < 10 AU, and for one that can decay faster than it 
would be disrupted by two-body relaxation it is more like 5 
AU. At a smaller pericentre distance than is given by this 
expression, gravitational radiation dominates the evolution; 
conversely, at a larger pericentre distance, two-body relax- 
ation dominates. 

At the MBH masses ~ 10 7 M Q that we consider, there 
may or may not be time for the stars to relax dynamically, 
hence it is not clear which value of a to take. If strong mass 



segregation occurs then a = 2 is likely (|Alexander & Hop- 


man||2009 


Preto & Amaro-Scoanc 2010| |Amaro-Seoane & 


Preto 


2011 


), but flatter slopes may also be relevant, par- 



ticularly if there has been scouring by a previous massive 
black hole merger and the system has not yet readjusted. 
For a selection of slopes we find 
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r p , crit « 14 AUM 7 13/10 (a/1000 AU)" 3/5 , a = 3/2 

r p ,crit « 7 AUM 7 5/4 (a/1000 AU) _1/2 , a = 7/4 

r p ,crit « 4 AUM 7 6/5 (a/1000 AU)" 2/5 , a = 2. (6) 

We will simplify by assuming that gravitational radia- 
tion is unimportant until r p — r p , cr i t , at which point it takes 
over completely with no further influence from two-body ef- 
fects. If this is true, then the next question is whether r p , cr it 
is greater than the tidal radius. If we focus on main sequence 
stars of mass m < Mq, then over a wide range of masses 
their radii are reasonably fit by i?* « 0. 85 Rq^/Mq) 2 ^ 3 
( Demircan Sz Kahraman|1991 I and the tidal radius is 



R 



1/3 



1.3 AUA'/ 7 1/3 mJ /3 



(7) 



Thus we see that stars in this mass range will typically enter 
the gravitational radiation regime before they are tidally dis- 
rupted. Given that the critical pericenter is just a few times 
the tidal radius, and that many aspects of this calculation 
are uncertain, it is quite possible that although tidal effects 
drop off very sharply with distance they could have an im- 
pact on the orbit outside pr- An exploration of this possibil- 
ity would require careful hydrodynamic simulations, but for 
our purposes we will assume that they are not dominant. 

Assuming that this is the case, we can compute the 
eccentricity of the orbit at the point that the pericentre dis- 
tance equals pr, when (as we show in the next section) the 
star is likely to be tidally disrupted instead of settling into a 
phase of steady accretion. We calculate the eccentricity by 
noting that to lowest (quadrupolar) order, pure evolution 
via gravitational radiation conserves the quantity 



C = ae 



-12/19 



(1 



-, 121 2 

1 H e 

304 



-870/2299 



(8) 



( |Peters|[l964| ). We saw in § [5] that the initial eccentricity 
after tidal separation is nearly unity, so 1 + e ~ 2. From 
our assumptions we also know that r p = a(l — e) = r PjC rit. 
Finally, if we assume that at the tidal radius the eccentricity 
is er <C 1, so that ar ~ Pr, we get 



-12/19 ^ -| o 

a T e T « 1.8r p ,crit 



e T ~ 0.4 



-19/12 



(9) 



For our three slopes the eccentricity at the tidal radius is 
thus 



e T » 0.01 M 7 - bai/3au m; J/3a (a/1000 AU) la/2U , a = 3/2 

e T « 0.03 M 7 - 209/144 mJ 9/36 (a/1000 AU) 19/24 , a = 7/4 

e T « 0.07 M 7 - 247/180 mJ 9/36 (a/1000 AU) 19/3 °, a = 2 

(10) 

We now explore the consequences of the star sinking 
inside the tidal radius with this eccentricity, and argue that 
tidal disruption is the most likely outcome if the SMBH 
has sufficiently low mass. We then demonstrate that tidal 



disruption with a small eccentricity leads to a different light 
curve than the more commonly considered tidal disruption 
of a star on a parabolic orbit. 



4 HYDRODYNAMICS NEAR AND INSIDE 
THE TIDAL RADIUS 

Suppose that the star sinks gradually under the influence of 
gravitational radiation towards the tidal radius. The tidal 
stresses increase as ~ (i?*/r) 6 , where R* is the stellar ra- 
dius and r is the distance from the SMBH. Therefore the 
star will be flexed and distorted, and internal modes will 
be excited as it sinks (for a recent discussion and simple 
model of this complicated process, see Ogilvic 2009). If the 
energy from these modes could be dissipated then the or- 
bit would undergo tidal circularisation and might end up 
in a stable mass transfer state. However, the energy that 
must be dissipated is significantly larger than the binding 
energy of the star. To see this, note that at the tidal radius 
Pr, we have tt = (3M/m) 1//3 R+. The binding energy of the 
star is w Gm 2 /R*. The binding energy of the orbit is 
Borb ~ GMm/rr- Circularisation of an orbit with eccen- 
tricity e at constant angular momentum releases an energy 
e 2 £ ; orb, so the ratio of released energy to stellar binding en- 
ergy is 



E, 



-1/3 2 



2/3 



(11) 



If M ~ 10 m the ratio is therefore > 



:10 4 e . From the 

previous section we found e ~ 0.01 — 0.07 for M = 10 7 Mq, 
so the energy required to circularise the orbit would be ~ 
3 — 150 times the binding energy of the star. If this energy 
could be released slowly this would cause no problems (note 
for comparison that in its lifetime the Sun will radiate a few 
hundred times its binding energy). However, the thermal 
(Kelvin-Helmholtz) time for solar-type stars is a few tens of 
millions of years, much longer than the inspiral time in our 
case and thus the tidal stresses will build up more rapidly 
than their mode energy can be radiated. 

The competition is therefore between the time needed 
for gravitational radiation to move the star into the tidal 
radius (where mass transfer will ensue) and the time needed 
for circularisation due to tidal dissipation to deposit a stellar 
binding energy into the star and thus, presumably, to tidally 
disrupt the star. Note that Alexander & Morris (20031 dis- 



cussed how tidal energy could produce "squeezars" with a 
different appearance from normal stars, without destroying 
the stars if the pericentre distance is sufficiently large. Here 
we are interested in the conditions for tidal destruction. 
To evaluate this we adapt the expressions from Lecontc 
(20101 for the energy deposition rate of tidal dissi- 



et al 



pation in a planet due to its eccentric orbit around a star. 
They find 



where 



-Etid 



N(e) = 



2Ki 



N a (e) 



N\e) 



1 + 



15 2 
2 c 



+ fe 4 + 



(1 



(12) 



(13) 
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N a (e) = 



1 + £ e 2 + ?f e 4 + We 6 + He 8 



(1 



16_ 

. e 2US/2 



l+§e 2 + §e 4 



(1 



(14) 



(15) 



and 



_i /Gm 2 



M\ 2 /7?„ N " 



1/2 



In the last equation Q is the quality factor of the star, 
a standard parameterisation of the rate of tidal effects on 
to the star. The magnitude of Q is notoriously unc ertain; 
values of Q = 10 5-6 are commonly used (see, e.g., Miller 
et al.||2009 for a recent example). If we use the expression 
(Zt = (3M/m) 1/ ' 3 i?* for the tidal radius, this last expression 
reduces to 



Gm' 

~r7 



(GM_\ 
\ « 3 J 



1/2 



(17) 



In the limit e <C 1 we find N a {e) ' 

Thus 



1 + 27e 2 , and £2(e 

-Etide 



: 1 + 
(Gm 2 



15 2 
2 c 



7 „-i 2 

4 Q 6 V 



1 + 23e 2 , 7V 2 (e) 

m"" 



(18) 



Thus the time needed to circularise the available energy ~ 
e 2 (M/3m) 2/3 (Gm 2 /-R.) at a = a T is 

T circ , tidc ^e 2 (M/3m) 1/3 (Gm 2 /i?„) 

x[(7/4)Q- 1 e 2 (Gm 2 /R t )(GM/a 3 ) 1/2 }- 1 
= 3 x 10 7 s QM 2/3 m~ 1/6 

(19) 

where in the second line we have substituted R, = 
0.85-Rom 2 / 3 and ht = (3M/m)R r . The circularisation time 
from gravitational radiation alone, at e <C 1, is 



,gw w (15/304)c 5 a 4 /(G 3 [iM 2 ] 
« 6 x 10 11 s 2/3 mJ /3 



(20) 



( Peters|1964| ), where in the last line we again substituted in 



rc, Gw/Tcirc,tide 



2 x 10 4 Q^ 1 Af~ 4/3 m 1 



4/3 1/2 







a — ar- Thus T c ir 

which for Q ~ 10 5_b is typically less than unity, hence only 
a fraction T c i r c,Gw/Tcirc,tide of the circularisation energy will 
go into tidal heating. Note, however, that for lower masses 
the eccentricity at the tidal radius is larger (scaling roughly 
as M -3 / 2 for our three power laws) and that the ratio of 
circularisation energy to the internal binding energy scales 
as e 2 , meaning that the total energy dissipated tidally scales 
as ~ M~ 4 , approximately. Thus even for Q = 10 6 , several 
times the stellar binding energy will be dissipated for M < 
3 x 10 6 M Q . 

If instead the SMBH mass is large, so that gravitational 
wave circularisation dominates over tidal circularisation, we 
expect that the star will settle into a period of steady mass 
transfer. The rate would be such that it balances the inward 
movement due to gravitational radiation, i.e., the charac- 
teristic time would be of order the gravitational radiation 
time. For our typical values, this is roughly 10 years, im- 
plying a rate of ~ 10 -5 M@ per year. Even if the luminosity 
is produced with an efficiency of 10%, this would produce a 
luminosity of only ~ 10 41 erg s _1 , weak enough and steady 
enough that it would not be distinguishable from a standard 
low-luminosity AGN. We therefore focus on the possibility 



that the star is tidally disrupted and that its debris is sub- 
sequently accreted by the SMBH. 

If a star is disrupted from a low-eccentricity orbit the 
evolution of its tidal debris proceeds differently than if it is 
disrupted from a parabolic orbit. To see this, note that in 



the original argument of Rees ( 1988 1 it was demonstrated 



that the spread in the binding energy of the debris is com- 
parable to the range in orbital binding energy from one side 
of the star to the other. If m/M ~ 10 -7 , therefore, the frac- 
tional spread is ~ (m/M) 1/3 ~ 10~ 2 . As a result, if an orbit 
with a pericentre r ~ rr has an eccentricity e > 0.99, the 
debris semi-uniformly samples binding energies from zero 
to the binding energy of the original stellar center of mass. 
The assumption of exactly uniform sampling (equal mass for 
equal range in binding energy) leads to a mass accretion rate 
that scales with time t as t~ 5 ^ s ; this law is more generally 
obtained at late times even for more realistic assumptions 
about stellar structure (e.g. Lodato et al.|2009 l. In contrast, 
if the spread in debris energies is much less than the average 
binding energy (corresponding to c < 0.99 in our exam- 
ple), then to lowest order the debris moves in a thin stream 
that intersects itself and settles within a few orbits into a 
remnant disc. 

Jl990|, who 



Such discs were studied by Cannizzo et al 



found that for plausible opacities the accretion rate would 
decay more gradually, e.g., M oc t -1 ' 2 for Thomson scatter- 
ing. Moreover, because the debris would all be bound to the 
SMBH (unlike for the parabolic case, where roughly half the 
stellar mass escapes to infinity), the accretion rate could be 
quite substantial for comparatively low-mass SMBHs. For 
Thomson scattering, the expressions from |Cannizzo et al.| 
(1990b lead to 



M = 2 x 10 23 g 



-1 ^_£^ 4 / 3 _7/9 



5/3 



(21) 

Here p is the average density of the star in units of 
g cm -3 , AM is the mass of the remnant disc (which will 



initially be the mass of the star) and a is the Shakura & 
Sunyaev ( 1973 1 viscosity parameter. For Mj < 1 this there- 



fore has the possibility of shining at luminosities that are 
a significant fraction of the Eddington luminosity Le = 
1.3 x 10 45 M 7 erg s" 1 assuming an efficiency L/Mc 2 = 0.1. 

As pointed out to us by E. S. Phinney (2010, personal 
communication), depending on the very uncertain details of 
how tidal energy is deposited, is it possible that there will 
be a gravitational wave signature that attends the electro- 
magnetic signature of disruption. In particular, it is not well 
established whether the tidal energy is deposited uniformly 
in the volume of the star or primarily where most of the 
matter is (both of which would lead to full disruption) or 
primarily in the envelope. If the last occurs, then the enve- 
lope would be stripped and lead to significant accretion with 
the characteristic decay discussed above, but the dense core 
would survive and could spiral in further. This would lead 
to a coincident gravitational wave signal that could be de- 
tected with the proposed LISA if the source is close enough 
( Freitag||2003[ ). 
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5 CONCLUSIONS 



In his work, Hopman ( 2009 1 estimates that for a galac- 



tic nucleus such as ours, the tidal separation rate of bi- 
naries which start far away from the MBH is 



pGC 

tid sep 



7 x 10 _7 (/b| uc /0.05)yr~ i , where "GC" stands for Galac- 
tic Center and /;, is the fraction of stars in binaries. Fig. (6) 
of |Hopman | ( |2009 |l shows that the rate increases when we go 
to higher energies, because the loss-cone is depleted, allow- 
ing more binaries to ''survive" in their way to the GC. Yu fc 
Tremaine ( 2003 1 estimate that the number is enhanced by 



an order of magnitude by binaries not bound to the MBH. 
More remarkably, the event rates can be at least temporarily 
enhanced by many orders of magnitude if one considers the 
role of massive perturbers, such as giant molecular clouds 
or intermediate-mass black holes, which can accelerate re- 
laxation by orders of magnitude as compared to two-body 
stellar relaxation (Perets et al. 20071. Another important 



potential boosting effect is the possibility that the potential 



is triaxial and not spherically symmetric (Poon & Merritt 
20021 120041 iMerritt & Poon|2004|). Taking these effects into 



account, we assume r^ sep ~ 10 _s (/ b | /0.05) yr" 1 . The 
fraction of main sequence stars that will eventually spiral 
into the SMBH after tidal separation is at least a few per- 
cent, so a plausible estimate of the total event rate for tidal 
disruptions of a single star originated by a separated bi- 
nary in a Hubble time is r GC ~ 10~ 7 (/ b | GC /0.05) yr" 1 , and 
it could be higher. This rate is probably a subset of the 
rate at which single stars are likely to encounter SMBHs 
on parabolic orbits (see Amaro-Seoane et al.|2007 for a dis- 
cussion of such extreme mass ratio inspirals). It is therefore 
possible that events with the L oc t~ 1,2 decay characteristic 
of low-eccentricity disruption may have rates smaller or sim- 
ilar to events with the L oc t~ 5 ^ 3 decay that is expected to 
be signatures of disruption of single stars in galactic nuclei 
and that is consistent with the initial decay of the recent 
Swift event Sw 1644+57 (Bloom et aLpmj ) . 
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